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Tliis  report  describes  the  capabilities  and  operation 
of  CONIXXIR  II,  a revision  of  the  surface  fitting  and  mapping 
prognim  originally  developed  by  Mr.  A.  Beharcll,  now  at  the 
University  of  t’algan'.  The  program  is  designed  to  fit  a 
geometric  surface  in  three-dimensional  space  to  any  given 
set  of  coordinate  data  points.  The  computed  "smooth"  surface 
is  then  displayed  by  drawing  specified  contours  on  the 
Talcorn]-)  93(i  driuii  plotter.  Although  this  nvinual  includes 
material  contained  in  the  University  of  Talgar)'  report  on 
UON'RXIR,  emphasis  lias  been  placed  on  the  input  changes  and 
additions  made  to  facilitate  tlie  use  of  this  [irogram  on  the 
nxi  bOnO  computer  systems  for  engineering  calculations  of  fluid 
d\oi;imics  [iroblems.  Details  arc  presented  of  the  control  cards 
ami  input  data  preparation  for  operation  of  the  prognim.  A 
sam^ile  problem  is  included. 

I.  INTROniO'lON 

Hie  contouring  teclinique  to  be  described  was  developed  by 
Mr.  A.  Beharell  and  later  modified  by  him  at  the  University  of  lalgan'. 
llis  program  was  purchased  by  the  Comiuiter  facilities  Division  of  tlie 
Uomjiut at  i on  and  Matliematics  IVpartment.  Because  of  extensi\’e  interest  in 
a general  putqiose,  versatile  contour  mapping  program  in  the  Numerical 
fluid  'lechanics  Branch  (Code  1843),  the  tape  containing  the  code  .md  a 
cop>'  oi  the  User's  M;inual  were  obtaiiieil.  Ilie  progr;im  was  then  modified 
to  more  efficiently  meet  the  needs  of  several  projects  in  the  Branch 
which  re(|uired  contour  diagriims.  lor  e.\;imjile,  the  IXIVIVUIR  11  jirogr.u-!  wa-- 

* "UON'DXIR,  A Surface  fitting  and  Mapping  Program,"  fhe  Universitv  of 
(Xilgary,  Itepartment  of  Computer  Sciences,  fSUCi-l'JD  (febniarv  19"bl. 
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used  to  obtain  pictures  of  wave  patterns  generated  by  an  Air  Cushion 
Vehicle,^  an  example  of  which  is  shown  in  Figure  1.  Another  application 
of  this  versatile  plotting  program  has  been  made  in  the  area  of 
structural  stress  analysis. 


Figure  1.  Waves  Cicnerated  by  an  Air-Cushion  Vehicle 


Haussling,  Henry  .1.  and  Vanliselt ine,  Richard  T. , "Unsteady  Air-Cushion 
Vehicle  Hydrodynamics  Using  Fourier  Scries,"  d.  of  Ship  Research,  Vol.  20 
No.  2,  .June  1976,  pp.  79-84, 
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A major  chiuige  made  in  the  original  CONTOUR  program  involved 
dividing  the  program  into  "tasks"*  to  reduce  the  computer  core  require- 
ments. In  addition,  the  fixed  input  arrangement  was  changed  to  free 
fonnat  for  ease  of  use.  An  extra  input  card  was  added  to  provide  several 
additional  options. 

Hie  resulting  program,  which  we  call  CONTOUR  II  to  distinguish  it 
from  the  original  version,  is  a surface  fitting  and  mapping  program  which 
complements  a number  of  other  graphics  tools  developed  in  the  Branch. 

I'hcse  other  programs  make  use  of  both  passive  and  interactive  graphics  i 

and  combinations  of  the  two.  for  example,  X'YZPLOT'^  is  a passive  graphics  i 

program  for  displaying  the  results  of  potential  flow  calculations.  : 

4 = 

HNCrhOr  is  a computer  program  subroutine  for  generating  graphs  of  data  \ 

in  fonnats  wtiich  are  useful  for  engineering  studies.  DI-SICN^  was 

developed  as  an  interactive  computer  graphics  package  for  displaying  and  . 

modifying  projections  of  three-dimensional  representations  of  vehicle 

models.  III.SICN  has  since  been  modified  and  enhanced  into  a comprehensive, 

modular  progr:im  called  IMACI'.  IMVili  is  an  interactive  program  incor-  , 

porating  the  three-dimensional  features  of  DliSlGN  ,'ind  also  providing  two- 

dimensional  graphing  capability  for  representing  FORTRVN  arrays  of  data 

as  cross  plots.  St!  4()(i0  and  Calcomp  hard  copies  are  easily  obtained. 


* 

A task  is  a section  of  code  which  may  be  executed  independently. 

Morawski,  Paul,  "XYZPLOT;  A n-iree-lUmensional  Graphics  Package  for 
Muid  IKTiiimics  (Z'llculat  ions  ,"  ('om]xitat  ion  and  Mathenvitics  llepartment , 
Departmental  Report  OID- l.s-"'.S  (August  IP'-'i). 
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Marquardt,  Mar>'  Beth,  "PNGPI.OT:  ,\n  l.ngineering  Plotting  Progr;im," 
(Y)mputation  and  Miithematics  Department,  Technical  Note  CiMIi-9-"4 
(I'ebniary  19'4). 

Kelly,  Barbara  M.  and  M;irquardt,  M.iry  Beth,  "Interactive  Helicopter 
Ik'sign:  Geometr>’  Package  User's  Ntinail,"  liomj-Hitat  ion  and  Mithematics 
Department,  |)epartmental  Report  Q1D-Z8-74  (September  19~4). 

Haas,  Melvin  I..,  "Interactive  IMAGT.  Display  .System,”  I'orninitat  ion  ;uid 
Miithenvit  ics  Ik'partment  Report  (in  jirep.-irat  ionl  . 


Any  series  of  data  points  which  can  be  represented  by  X,  Y,  and  Z 
coordinates  (or  digital  triplets)  can  be  used  to  describe  a surface  in 
three  dimensions.  A contour  is  the  line  connecting  points  of  equal  value. 
The  CONTOUR  II  program  is  designed  to  fit  a surface  to  any  series  of 
digital  points  (X,  Y,  Z) . This  surface  is  then  displayed  for  visual 
examination  by  drawing  the  specified  contours. 

The  surface  is  generated  by  first  calculating  a series  of  planes 
through  the  data  points  using  various  weighting  factors  and  then  inposing 
a grid  over  the  planes.  It  is  assumed  that  the  result  of  this  operation 
is  a continuous  surface.  The  contours  are  then  smoothed  before  plotting 
to  provide  a more  realistic  picture. 

This  type  of  program  has  been  used  in  many  areas  where  graphical 
expression  is  helpful  in  data  analysis.  Some  examples: 

• Meteorology  - barometric  pressure  and  other  climatic  factors. 

• Geology  - topographic  and  subsurface  mapping. 

• Civil  Engineering  - location  of  structures  with  respect  to 

such  areas  as  topography. 

• Oceanography  - water  salinity. 

Space  Exploration  - representation  of  magnetic  field  intensity. 

• X-Ray  Crystallography  - electron  density. 

• Medical  X-Ray  Technology  - radiation  intensity  within  the 

human  body  as  a result  of  x-ray  treatment. 

■ Air  Pollution  - pollution  indices. 

• Fisheries  - catch  of  different  types  of  fish  within  a given  , 

time  interval. 

• Census  - census  information,  population  and  related  itenrs. 

Fluid  Mechanics  - surface  elevations,  equi -potential  lines, 

streamlines,  pressure  distributions. 

Much  of  this  report  is  taken  directly  from  the  University  of 
Calgary  report.  Changes  have  been  made  mainly  in  the  input  description. 
Control  Card  examples  reflect  the  ITI’NSRIXl  CIX'  6000  operating  system.  The 
original  Calgary  version  was  inoperable  on  the  computers  at  our 
installation. 
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1 1 . SYSTIW  DILSCRI PT  ION 
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I nn-;  Minnon  of  coNiaiRiNG 

I In  order  to  take  advantage  of  the  flexibility  of  control  over  para- 

\ meters  involved  in  the  contouring  process,  the  user  should  be  familiar 

i with  the  method  of  transforming  data  into  a contour  map.  A contour  map  is 

a two-dimensional  representation  of  a three-dimensional  surface,  derived  by 
drawing  the  lines  of  intersection  of  a series  of  parallel  planes  with  the 
j three-dimensional  surface  (see  Figures  2 and  3). 

! In  the  description  following,  the  Cartesian  coordinate  system  will  be 

used.  It  will  be  assumed  that  the  scries  of  parallel  planes  is  parallel 
i to  the  X,Y-plane,  cmd  intersect  the  2-axis  at  regular  internals  as  sliown 

in  Figure  2. 


l igiiri'  2.  Sample  Ilirci'  Dimensional  Snriaee  ( lit  h\  Parallil  I’lines 


.3 
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Figure  .V  CuiKoiir  \l;ip  ProdtRinf  hy  VFippiug 
Inlerseetion  Lines  of  Pbnes  with  Siirtaee 


The  contouring  process  uses  a set  of  user  supplied  DATA  1X11  NTS. 
rhese  are  a series  of  points,  described  by  X,  Y,  and  Z coordinates  which 
define  the  three-dimensional  surface  to  be  mapped.  (Sec  Figure  4). 


Figure  4.  Data  Points  in  Three  Dimensions 


The  problem  is  to  create  a CONTOUR  SlIRFACH,  that  is,  a three- 
dimensional  surface  which  best  passes  through  the  given  DATA  POINTS,  ;md 
best  predicts  trends  between  the  data  points.  Approximation  of  the 
CONTOUR  SURFACF,  begins  with  the  determination  of  a PIANl'.  surface  at  each 


{) 


HATA  roiNT,  whidi  best  reflects  trends  suggested  by  the  closest 
neighboring  DATA  POINTS.  To  determine  the  PIANl:  surface  at  a particular 
DATA  roiNT,  consider  the  two-dimensional  map  of  the  DATA  POINTS  in  the 
X,V  plane  (Figure  5).  Tlic  area  around  each  DATA  POINT  is  divided  into 
i-XT.ANTS,  that  is,  8 equal  pie-shaped  pieces  (Figure  6).  In  each  of  these 
(Xli'ANTS  a search  is  conducted  for  the  closest  NF.IQIBOR  in  the  X,Y  phine 
to  the  DATA  POINT  being  considered  (Figure  7).  The  PIAN'i;  surface  is 
chosen  to  pass  through  the  DATA  RMNT  ;ind  to  pass  "most  closely"  through 
the  NFIdlBORS,  with  a weighting  scheme  used,  such  that  the  PLANT  surface 
fits  more  closely  to  nearer  chosen  N'TIGIIBORS  (Figure  8).  Thus,  a PIANT 
surface  is  associated  with  each  given  DATA  POINT. 

Ay 
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Fiuiiri'  5 l):il.i  I’uinis  Mappiil  in  Iwu  Dimfiisiuns 


I lull u’  ^ Si'li'i  lmn  III  NcivliBiii s 
.ihiuti  .1  ( ■ i\  ni  I ).ii.i  I’ll  nil 


I iunri'  (>  OilanN  ahiini  a 
( liu'ii  Data  I’liint 


Figure  8.  Dclemiinalion  of  Plane  Surface  Using  Cliosen  Neighbors 


The  next  stage  in  computing  the  CONIOUR  SURFACH  superimposes  a 
regularly  spaced  grid  on  the  X,Y  plane  (Figure  9).  Fach  intersection  of 
the  grid  lines  defines  a GRID  POINT  at  which  a GRID  VALIIF  or  Z -value  is 
computed.  The  resulting  set  of  GRID  VALUES  is  assumed  to  define  the 
CONTOUR  SURFACE. 

The  Z -value  at  each  GRID  TOINT  is  computed  as  follows.  At  each  GRID 
roiNT  a search  is  made  for  the  N DATA  TOINTS  (where  \ is  a number  from  1 
to  20)  which  are  nearest  in  the  X,Y  plane  to  the  GRID  POINT  (Figure  10). 
The  PlANi;  surface  computed  earlier  for  each  of  these  DATA  POINTS  inter- 
sects a line  through  the  GRID  POINT  and  parallel  to  the  Z-axis,  giving  a 
Z-value  (Figure  11).  Tlie  N Z-values  thus  found  at  the  GRID  TOINT  are  each 
weighted  according  to  the  distance  of  the  DATA  fXUNT  from  the  GRID  POEM 
such  that  the  closest  has  the  greatest  weight.  Hie  weighted  average  of 
the  N values  of  Z is  taken  as  the  value  of  the  CONTOUR  SURFACE  at  tlie 
GRID  POINT.  ITiis  process  is  repeated  to  compute  a Z-value  for  eacli  GRID 
WINT  on  the  X,Y  plane. 


I'iuiiri'  Sii)K’rim)H)soil  (iriil 


l 10  Si‘k‘1  lion  lit  Neighbo 


I iuiiri'  I I Di'U'i iniiKil inn  nl  ( miltnir  V':iliio  ;il  :i  (iriil  Point 
IkiM'd  on  llu"  PLinos  ol  ilii'  Ni’iglihors 


Having  thus  dctc'mnncd  regularly  spaeed  values  of  tlie  fONTiXIR 
SlIRlACd,  over  tlie  entire  map  region,  the  next  step  is  to  represent  the 
three-dinK’nsional  surface  so  defined  in  two  dimensions  as  a contour  ruii'. 
Placement  of  tlie  contour  lines  on  tlie  ma]i  pioceeds  with  tlie  cons ulerat  ion 
of  the  ('OVliXlR  JJATLS,  one  at  a rime,  from  Iriwest  to  higliest  . for  a 


particular  CCWTOUR  LEVEL,  the  GRID  VALUES  determined  above  are  searched 
until  two  adjacent  GRID  POINTS  are  found  which  have  values  above  and  below 
the  particular  COIWIUR  LEVEL.  These  two  GRID  POINTS  define  a GRID  SQUARl: 
in  which  a contour  line  starts.  The  starting  point  is  determined  by 
linear  proportion  between  the  two  GRID  POIMTS  (Figure  12a) . This  GRID 
SQUARE  is  then  searched  for  another  side  whose  end  points  lie  above  and 
below  the  CONTOUR  LEVEL  under  consideration.  The  contour  line  is  made  to 
exit  from  the  current  GRID  SQUARE  to  an  adjacent  one  through  this  side, 
again  by  linear  proportion  between  the  end  points  (Figure  12b). 


Figure  12.  Determination  of  a Contour  Line  through  a Grid  S<|iLire 


The  particular  CONTOUR  LEVEL  line  proceeds  from  adjacent  GRID  SQU.-\RI'. 
to  adjacent  GRID  SC^ARE  until  it  closes  on  itself  or  reaches  the  edge  of 
the  map  region.  At  this  point,  another  starting  point  is  found  for  the 
CONTOUR  LEVEL  or,  if  these  have  all  been  exiiausted,  the  next  level  is 
begun.  In  order  to  improve  the  appearance  of  the  resulting  contour  lines 
a smoothing  function  is  applied,  so  that  a scries  of  straight  lines  and 
angular  comers  is  avoided. 
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DILSHA'INC;  A CONTOim  HM’ 

Several  factors  should  he  taken  into  consideration  before  attem}itin^ 
to  fill  out  a set  of  control  cards  for  the  contouring  of  a set  of  data: 

1)  In  order  to  introduce  a set  of  DATA  IMINTS  into  the  contourinj^ 
prognun,  they  must  be  assij^ned  X,  V,  iind  Z coordinates  using  some 
convenient  system  of  units  (DATA  UNITS).  ITie  maximum  number  of  DATA 
IXU.VrS  in  one  iruip  is  998. 

Z)  Hie  tmip  boundaries  iraist  be  ascertained  either  by  studying  the  rmaximum 
;uul  minimum  X ;ind  V coordinates,  or  by  exrutuning  an  existing  map  of  the 
area  which  shows  tlic  positions  of  the  DATA  TOIN'IS  (figure  1.^). 

-A)  A value  must  be  assigned  to  the  lower  left  corner  of  the  map  in  X and 
V DATA  UNITS,  fniis  is  the  X ORIUIN,  V ORIGIN'  on  the  control  card). 

4)  rhe  distribution  of  the  DA'l'A  POINTS  sliould  lie  studied  ;ind  the  size  of 
t)ie  GRID  S(^IARTS,  which  the  prognim  will  suiter inpose  on  the  nvip  surface, 
sliould  tic  determined.  A good  rule  is  to  tr>’  to  make  the  grid  size  such 
that  no  more  tlwm  one  DMA  IXUNT  falls  witliin  each  GRID  SQUARI  . However, 
this  is  not  always  practical,  since  the  data  may  be  denselv  concentrated 
in  some  areas  and  sparse  in  others.  In  a case  like  this  it  is  best  to 
choose  some  intermediate  value.  It  is  stronglv  recommended  that  square  or 
neaiT\'  scjuare  rectangular  cells  !ie  used.  Ihe  grid  size  is  defined  on  the- 
control  card  as  X GRID. 

.s)  Another  consideration  in  detennining  grid  size  is  the  tact  that  the 
program  will  only  liandle  bSOO  GRID  IMINTS.  Die  numlx'r  of  GRID  DUVIS  is 
determined  liy  dividing  the  total  nimiber  of  DMA  UMIG  along  the  X-axis  bv 
X GRID  VA!,Ui.  to  give  the  num)ier  of  GOfUMNS  I , ,ind  dividing  the  total 
number  of  data  units  along  tlie  l-axis  b\-  V tlRlD  \'\1.U1  to  gae  the  number 

of  ROKS-1.  rhe  [iraluct  of  the  ROWS  and  fOl.UMNS  gives  the  number  of 

GRID  I’OIN'IS.  In  general,  doulil  ing  the  nimiber  of  rows  aiwl  columns  will 
increase  computer  time  by  a factor  of  4 (1 igure  14). 

())  Ihe  GO.NTOIIR  INTTRVAT  may  be  detemuned  In  examining  tlie  range  and 
variation  of  tlie  Z coordinates  of  the  data.  Usually  a map  witli  denselv 
packed  contour  lines  is  not  pleasant  to  look  at  and  takes  much  longer  on 

the  coTiputer.  Also,  if  the  contour  interval  is  cliosen  too  fine  for  the 

vlata,  the  smoothing  algorithm  used  m.i\-  cause  some  contour  lines  to  crosc. 


■’I  llic  appearance  of  the  contour  map  produced  by  CONTCXJR  from  a given 
set  of  data  is  controlled  by  a number  of  parameters  in  the  contouring 
process.  These  parameters  can  be  modified  by  the  user,  as  desired,  to 
achieve  a satisfactor>'  nvip.  Control  of  CONTOUR  is  designed  so  that,  if  the 
user  does  not  select  a value  for  certain  parameters,  a standard  value  will 
be  assumed.  ITiis  set  of  standard  values  has  been  chosen  to  give  good 
results  with  most  sets  of  data.  The  user  should  also  note  that  contouring 
is  a process  of  approximation.  ITic  accuracy  of  this  approx invit ion  is 
affected  to  some  extent  by  the  choice  of  parameters  and  also,  of  course, 
by  the  accuracy  ;uul  density  of  the  data  points  supplied.  Some  increase  in 
accuracy  can  be  achieved  by  choosing  other  than  the  st:indard  set  of 
parameters.  Tlus  will  usually  cause  the  program  to  run  longer  and  hence 
be  more  c.x]iensivc  to  use.  The  following  notes  indicate  the  parameters 
which  the  user  may  alter.  Ilie  standard  value  for  each  is  show7i  in 
brackets  after  the  n:uiie. 

MINQ11;R  Of  (XT.WIS  (5) 

Ibis  is  the  minimum  number  of  oct;ints  around  each  DATA  POINT 
which  must  contain  a neighbor  before  that  point  is  considered 
in  the  contouring  process.  If  the  point  meets  this  requirement, 
a PIAM;  surface  will  be  fitted  through  it  which  will  eventuallv 
influence  the  determination  of  the  GRID  V.-d.UbS. 

UTdOirilNC.  TACTOR  i'OR  Tin;  DAi'A  IXUNTS  (-2) 

As  mentioned  in  the  description  of  the  contouring  mctluxl,  a 
weighting  sclieme  is  used  which  causes  the  PIAM.  surface  to  fit 
more  closely  the  chosen  neighbors  to  a DATA  IXUNT.  lor  the 
standard  value  of  -2,  this  weighting  scheme  is  roughly  equivalent 
to  using  the  inverse  square  of  the  distance  of  neighbors  from 
tlie  point  being  evaluatcil.  A value  of  -1  is  roughh'  equivalent 
to  linear  weighting.  If  a value  of  DD.D9!'  is  enteied  here,  the 
weighting  will  be  comidetely  eliminated.  Positive  values  and 
zero  are  not  recommended,  llie  effect  of  increasing  the  weight 
(making  it  more  negative)  is  to  force  contour  lines  to  more 
accurately  honor  the  giien  DAI  \ I’lMVls;.  However,  in  areas  ot 
s[iarse  data  an  increase  in  weight  otten  causes  the  appearance 


r 


of  extreme  highs  and  lows. 

NUMBER  OF  NlilGHBORS  UO) 

This  is  the  number  of  closest  neighbors  to  a grid  point,  which 
will  be  used  in  determining  a value  at  that  grid  point. 

WEIGKTING  FACTOR  FOR  THE  GRID  (-2) 

Used  in  a similar  manner  as  the  WEiaiTING  FACTOR  I'OR  TUT  DATA 
POINTS.  This  weight  is  applied  to  the  Z-values  calculated  for 
the  grid  point  based  on  the  PLANES  from  each  of  the  chosen 
neighbors. 

RANGE  (5.) 

This  is  the  radius,  in  inches,  within  which  the  closest  NUMliER 
OF  NEIGHBORS  must  be  found,  if  a grid  point  is  to  be  evaluated. 

8)  A number  of  parameters  are  also  available  to  control  the  plotted 
output  on  the  contour  map.  Some  of  these  parameters  may  be  omitted  in 
which  case  standard  values  will  be  assumed  by  the  program.  The  following 
paragraphs  suggest  some  of  output  variations  possible.  Conplete  details 
of  these  parameters,  including  the  standard  values  assigned  by  the  program, 
are  described  in  section  3. 

The  user  'nay  choose  to  include  or  omit  the  plotting  of  a title,  the 
axes,  a symbol  marking  the  exact  location  of  each  data  point,  the  2-valuc 
of  the  data  point  (which  may  be  placed  in  a variety  of  positions 
relative  to  the  exact  location),  and  the  values  of  the  contours.  He  may 
also  choose  to  have  some  contour  lines  drawn  heavier  than  others  and  may 
choose  to  omit  contours  of  less  than  a specified  length  or  based  on  less 
than  a specified  number  of  data  points.  Unless  the  user  overrides  the 
standard  choice  of  minimum  and  maximum  contour  lines,  by  entering  these 
on  the  control  card,  the  program  will  choose  these  limits  based  on  the 
maximum  and  minimum  Z-values  of  the  DATA  POINTS  (or  if  the  data  is  GRID 
POINTS,  the  maximum  and  minimum  GRID  POINT  values). 

A special  feature  allows  the  plotting  of  a single  set  of  contour 
data  as  two  separate  contou*'  maps.  This  is  called  MAI’  SF.P(\R\T  1 ON . ITie 
most  common  application  of  this  feature  would  be  in  the  case  of  rcsidiuil 
maps,  which  involve  Z-valucs  above  and  below  zero.  If  MAP  SliPARATION 


14 


were  specified  at  contour  level  zero,  two  contour  mxips  would  be  drawn. 
iTie  first  would  show  all  contour  levels  below  zero,  and  the  zero  line. 

Die  second  would  .show  the  zero  line,  and  those  above  zero. 

.Some  additional  options  j>iv'c  the  user  .several  useful  plotting 
capabilities,  (inc  of  these  features  is  referred  to  as  M\P  ,SUri:RIM}’0.sn'I(>\, 
which  allows  two  different  sets  of  data  to  he  plotted  on  the  s.ame  contour 
map.  Hris  is  useful  for  comparison  of  contours  or  for  elimination  of  a 
specific  contour  level.  I'or  cxamjrle,  if  it  was  beneficial  to  eliminate 
the  zero  contour  level,  the  u.ser  would  specify  negative  lower  ;ind  upper 
contour  limits  on  the  data  cards  to  produce  the  first  map  ;ind  jiositivc 
lower  and  upper  limits  on  the  nc.vt  set  of  data  cards  (with  the  MM' 
bl!Pi;i< IMPOSITION  flag  .set  to  1 ) . 

The  u.ser  may  also  specify  a "noise"  toleraiKC  based  on  data 
generation  errors.  With  tliis  option  the  user  c;in  eliminate  meaningless 
contours  from  the  map.  A default  is  used  if  no  value  is  specified;  tiiat 
means  a tolerance  value  is  preset  in  tlic  program  if  none  is  specified  in 
tlie  input. 

Ibe  user  may  request  dashed  lines  for  contours  representing  negative 
values  by  specifying  a dash  patteni  and  a dash  lengtii  in  inches.  A 
default  length  is  already  set  in  the  program. 

Body  outlines  may  be  dravoi  on  the  map  by  inputting  the  file  number 
wiierc  the  user  has  set  up  X,Y  points.  Oaps  or  "holes"  may  be  specified. 
Outline  points  must  he  in  user  data  units. 

ixnrr  oitiovs 

rive  t>pes  of  input  are  available  in  the  (UNTOUR  prognim.  Hiese 
arc  described  as  follows: 

T>pe  I Ibis  is  the  standard  t\pe  of  input  most  comnKmlv 

used:  Die  X.V.Z  coordinates.  Ibis  infoimiation  can 
be  read  from  am’  file  containing  one  set  of  values 
(ler  record,  (be  fonriat  of  the  data  mav  be  read  in 
free  fonn.it  or  as  specified  bv  the  user,  (be  end 
of  the  ini'ut  data  is  indicated  by  ;i  final  record 
cont. lining  the  di.-iracters  po. 


1' 


Type  2 


Type  3 


<! 


This  is  for  input  in  the  form  of  the  coefficients 
of  a surface  equation  of  order  1,  2,  3,  or  4.  Tlie 
number  of  coefficients  required  is  3,  b,  10  and  15, 
respectively,  corresponding  to  the  following 
formulae: 

First  Order 


Z C^+C^X+CjY  ^ 

Second  Order 

Z=C jH;  ,X+C^Y+C^X"+C^XY+C^Y^ 


Third  Order 


.3, 


..■) 


z=c^+c,x+c^y+c^x“h:^xt+c^y  +c-x'-k:^x  Y+c^xT-^c^yV 
Fourth  Order 

Z=C^+C2X+C,Y+C^X^+C^XYH:^Y^+C.X'^H:j^X^’+CqXY"+  Cj^^Y'"’ 

4 3 2 2 5 4 

+C, ,X  +C, ,X  Y+C, ,X  Y +c,  AT  +C, A 
11  IZ  1:>  14  Is 


Input  to  CONTOUR  consists  of  the  order  of  the  sur- 
face, number  of  coefficients,  and  a string  of 
coefficients  in  a free  format  (see  section  3 for 
format).  No  end  of  data  specification  is  reeiuired 

This  type  of  input  specifies  the  GRID  VALUliS 
directly.  Normally  these  will  have  been  generated 
in  a previous  CO.NTOUR  nin.  Input  data  consists  of 
the  number  of  rows  and  columns  of  the  grid  and  the 
set  of  values  in  free  format  (see  section  3). 

ilie  order  of  the  GRID  TOINT  values  is  row  by  row 
from  left  to  right  for  each  row.  Neither  input 
format  nor  end  of  data  indication  are  required.  A 
value  of  -1F..50  will  be  interpreted  the  s;une  as  a 
blank  value  and  no  contours  will  be  drawn  based  on 
such  a point.  The  grid  input  is  identical  to  tlie 
output  produced  by  the  PUNCH  grid  option.  Use  of 
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this  option  means  a subsequent  run  will  not  be 
required  to  recalculate  fpRlD  VALUliS. 

Tvpe  4 ITiis  option  uses  X,  Y,  Z data  and  the  phine 

coefficients  generated  from  them  in  a previous  run. 
lYiis  option  may  lx;  useful  when  it  is  desired  to 
produce  more  than  one  contour  map  for  a set  of  data. 
This  may  be  desirable  to  obtain  the  optional  set  of 
control  parameters.  Data  for  this  option  is 
produced  by  the  U7UT1;  CObFnCIHNTS  output  option. 
iTie  unit  number  of  the  file  must  be  specified  as  one 
of  8,  9,  or  10. 

T\pe  5 This  option  is  used  to  input  a set  of  GRID  POINTS 
generated  by  a previous  CO.NTOUR  run  using  the 
IvTUTF  GRID  output  option.  Input  par;imeters  required 
for  this  option  are  similar  to  T>pc  4.  It  may  bo 
useful  to  save  ;ind  restore  GRID  TOINTS  if  changes 
are  to  be  made  which  affect  the  plotting  of  tlic 
contours.  For  e.Xiimple,  the  second  nup  could  have  a 
different  contour  interval.  With  this  option, 
however,  it  is  not  possible  to  annotate  the  original 
data  points  as  they  are  not  preserved. 

(xrri’irr  options 

Several  output  options  are  available  for  operations  such  as  printing 
the  data,  pimching  cards, or  writing  to  a file.  Different  t\j)cs  of  data 
can  be  output.  Following  is  a description  of  the  available  options: 

PRI.NT  DATA  Any  t>T>e  of  input  c;m  be  printed  as  required. 

PRINT  COF.FFUT PINTS 

Ihis  lists  the  coefficients  of  the  PIANF  surfaces 
fitted  through  each  individual  DATA  !\''1NT.  Ihe 
data  is  printed  in  a sorted  sequence  of  N,V,Z. 

(Rata  IS  sorted  on  the  V coordinate.) 
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PRINT  GRID 


This  lists  the  entire  grid.  The  grid  numbers 
produced  begin  at  the  lower  left  comer  and 
continue  across  each  row  from  left  to  right. 

WRITE  COEFFICIENTS 

This  option  creates  a file  of  DATA  TOINTS  and  their 
associated  PLAN!:  coefficients.  /\n  identification 
record  is  also  written  at  the  beginning  of  the  file. 

Several  sets  of  COIiFFICILNTS  may  be  written  on  one 
file.  This  output  can  then  be  rewound  and  input  to 
succeeding  CONTOUR  operations  as  INPUT  TiPI.  4. 

WHITE  GRID 

This  option  creates  a file  of  GRID  POINTS.  Several 
sets  of  GRID  POINTS  may  be  written  on  the  same  file. 
This  file  becomes  INi’in'  ITPE  5 in  following  runs. 

PUNCH  GRID 

This  punches  GRID  TOINTS  in  the  format  described 
for  INPUT  TiPE  3.  This  feature  can  be  used  to 
modify  the  GRID  points  calculated  by  the  program. 


III.  CONTROL  CARD  DESCRIPTION'S 

The  control  cards  required  to  use  the  CONTOUR  program  consist  of 
three  types;  system  control  cards,  CONTOUR  parameter  cards,  ;ind  data  if 
required. 

SYSTIW  CONTROL  CARDS 

Following  are  two  examples  of  the  control  cards  required  to  nui 
CONTOUR  II  at  DTNSRDC  on  the  (TX’  6000  system.  For  each  run  the  user 
must  supply  his  own  tape  (slot  or  bin  number)  to  record  the  Calcomp 
plotting  instructions.  After  the  run,  the  user  must  submit  a request 
slip  to  have  the  tape  processed  for  Calcomp  hardcopy  output. 
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1)  n-iis  is  the  set  of  control  cards  required  if  X,Y,Z  coordinates  exist 
on  cards. 

CM .CMfiOOOn  ,T180  ,P2  .NTl'l . 

Cl  IMICI ; , xxxx , nnnnnnnnnn . 

ATTAQ 1 , CONTIIR , CWCWrOURTSK , I n=CAMV . 

VSN(TAPhll=CBxxxx')  bin  tape 

Ri;(piST,T;\Pi;ii,m,KiNC.. 

CONTUR. 

"89 

Projiram  parameter  cards 
X,Y,Z  data  cards 
PMIMM’ 

9 

(’"89 

11  is  the  lORTR'VN  unit  number  to  which  the  program  writes  the 
Calcomp  plotting  instructions. 

2)  Hiis  ex;imiilc  assumes  the  input  data  exists  on  a file  CWhCLiA'ATIONS 
;ind  is  to  be  read  by  CONTOini  on  FORIRW  logical  unit  9 (CIX:  file  TAPi.9). 
Hie  program  control  cards  (tNjies  1-b  including  ILMMAP  ;ind  last  card  (9)1 
are  on  a file  called  CONTROL. 

CM AICT, 0150000 ,i'180 , P2  .MLl . 

Cl  IMICI ; , xxxx , nnnnnnnniui . 

A'n'ACll  ,C0NT  .CLMRCONTODRl'SK , I IMLMR' . 

ATTACi  1 , 1APH9 .CAVliLLLVATIONS , I n=CAVL . 

ATTAG 1 , CONTROL  .CALMCONn'ROLCRDS , 1 n=CAHM. 
VSNrrAPL,ll=SL0T13=CALAlIlC.)  slot  tape 
RLQIILST,TAPL11,11I  ,RINC. 

C0N'r,C0NTR01.. 

”89 


I'ORTRW  units  8,  9,  and  10  may  be  used  for  input  as  idX'  files 
IAPL8,  TAI’LO,  and  TAPLIO.  Any  files  defined  as  direct  access  must  be 
AlTAClled. 

Hie  program  intenially  uses  files  LM'T.92  and  1 API  93.  Ihesc  are 
used  principally  for  saving  tcmporar\-  grid  files  .and  must  not  lie 
disturbed  bv  the  user. 


CONTOUR  PARAMETER  CARDS 


There  are  six  parameter  cards  required  by  CONTOUR.  These  are  Map 
Title  Card,  Contour  Map  Parameter  Card,  INPUT/OUTPUT  Card,  Processing 
Parameter  Card,  Plot  Parameters  Card  and  Extra  Options  Card.  All  cards 
must  be  included  in  the  specified  order.  All  cards  except  the  Map  Title 
Card  are  in  Free  Format. 

Free  format  (or  List  Directed)  input  data  consists  of  a string  of 
parameter  values  separated  by  one  or  more  blanks,  a comma,  or  a slash. 

When  the  value  separator  is  a slash,  remaining  list  elements  for  that 
card  are  treated  as  nulls,  that  is,  the  values  set  in  the  program  are 
used. 

To  repeat  a value,  an  integer  repeat  constant  is  followed  by  an 
asterisk  and  the  constant  to  be  repeated.  Variables  of  differing  t\q)es 
should  not  be  repeated  with  the  same  constant.  For  example,  to  read 
X,Y,I,J  the  user  must  input  2*7, 2*7,  and  not  4*7  (assuming  X and  Y are 
real  and  I and  J are  integer).  For  more  information  on  free  format  input, 
see  CONTROL  DATA  FORTRAN  Extended  Version  4 Reference  Manual. 

Map  Title  Card  This  card  specifies  a title  for  the  map.  Tlie  title 

is  printed  as  the  output  heading  and  plotted  above 
the  map . 

Col.  1:  1 

Col.  2-79:  Any  desired  title  information 

Contour  Map  Parameter  Card 

This  card  contains  variables  which  are  related  to 
the  size  and  scale  of  the  contour  map.  There  are 
no  defaults  for  any  field  on  this  card. 

Parameter  1 : 2 

Parameter  2:  X-ORIGIN  - the  value  of  X at  the  lower  left  in 

DATA  UNITS.’ 

Parameter  3:  X-GRID  - the  number  of  data  units  along  the  side  of 

a GRID  square  parallel  to  the  X axis. 

Par:imeter  4;  Y-ORIGI.N  - the  Y-value  in  data  units  at  the  lower 
left  comer  of  the  map. 
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Parameter  5: 


Parameter 

Parameter 

Panimeter 

Parameter 

Par:imeter 

iN'i'invoirrpirr  c: 

Par;uncter 

Parameter 


] 

Y-GRII)  V.-MAII;  - the  number  of  data  units  along  the  ■ 

side  of  a GRID  square  parallel  to  the  V-axis.  ' 

b:  CONTOUR  INTHRV.-M.  - the  number  of  Z data  units 

between  successive  contours. 

7;  COLUMNS  - the  number  of  columns  in  the  GRID  to  be 

superimposed  ov'er  the  map  area,  fhis  number  minus  1 
times  the  X-GRID  SIZL  determines  the  map  length  in  < 

inches.  | 

8;  X-GRID  SIZL  - the  length  of  the  X side  of  a GRID 

square  in  inches . ' 

9:  ROWS  - the  number  of  rows  in  the  GRID  to  be  super-  j 

imposed  over  the  map  area,  ibis  number  minus  1 ^ 

times  the  Y-GRID  SIZI:  will  give  the  height  of  the 
map  in  inches. 

10;  Y-GRID  SIZL  - the  number  of  inches  along  the  Y side 
of  a GRID  square. 

ard  Ibis  card  is  used  to  select  any  desired  injuit  ;md 
output  options  and  also  to  describe  the  input  data. 

1 ; 

Z;  iNPirr  iTPi: 

= 1:  Data  records  in  X,Y,Z  coordinates.  Records 
may  be  free  format  or  variable  foimiat  specified  as 
input,  l.nd  of  coordinate  values  is  indicated  by  a 
99  on  the  last  record. 

= 2;  Data  records  are  trend  surface  coefficients 
input  in  free  format. 

= ■>:  Data  records  arc  GRID  POINTS  in  free  fonaat. 

= 4;  Data  records  are  X,Y,Z  coordinates  and  pi. me 
coefficients  on  a file  generated  b\-  the  WRIII 
COl  IITCII'NTS  option  i i'ar;imetei'  1.7) 

= .S;  Dat.a  lecords  are  GRID  IXUN'l'S  on  a file  genei'.ited 
by  tlie  WKlli  GRID  option  (Par.imeter  Ml. 
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Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

F’arameter 

Parameter 

Panimeter 

f’arametcr 


3:  INl\rr  UNIT  - Fortran  logical  unit  of  the  input  data. 

The  card  reader  is  5,  Files  may  be  on  any  of  the 
units  8,  9 or  10. 

4;  RFWIND  - Enter  a 1 if  the  input  unit  is  to  be 

rewound  (positioned  at  start)  before  reading.  The 
card  reader  cannot  be  rewound. 

5:  INRff  JOB  - this  applies  only  to  Input  l^pes  4 and  5 

where  the  input  file  was  created  by  a previous 
contour  run  or  step.  It  is  the  number  of  the  job  in 
which  the  data  was  created.  It  is  a sequential 
number  corresponding  to  map  sequence.  This  number 
is  printed  at  run  time  for  each  job  step. 

6;  X TOSITION  - this  specifies  whether  the  X variable 
is  the  first  (1),  second  (2)  or  third  (3)  variable 
on  input.  The  default  value  is  1. 

7:  Y POSITION  - this  specifies  the  position  of  Y 

coordinate  as  1,  2 or  3.  Default  value  is  2. 

8;  2 POSITION  - similar  for  2-coordinate  position. 

Default  is  3. 

9:  ILND  n;ST  POSITION  - the  value  0 or  1 to  indicate 

that  the  position  of  end  of  record  characters  99 
and  read  with  an  A2  forniat  is  cither  the  first  (0) 
or  last  (1)  variable.  This  cntiy  reciuircJ  only  for 
Input  IN-pe  1, 

10;  PRINT  DATA 

= 1 to  list  the  input  data 

= 0 otherwise  I default) 

11 : PRINT  COEFFICII-NTS 

= 1 to  print  the  PIANT  coefficients 
= 0 othci'wise  (default) 

12:  PRINT  GRID 

= 1 to  list  the  GRID  IXMNTS 
= 0 otherwise  (default) 


Parameter  13: 


Parameter  14; 


Parameter  15: 


Parameter  16: 


WRITi;  COnr-r- U; I HNTS  - l-ortran  logical  unit  on  which 
coefficients  are  to  be  written  if  this  is  desired; 
valid  entries  are  8,  9 or  10.  fDefault=0,  no 
coefficients  written.) 

bTinT,  (.iRin  - F-ortran  logical  unit  for  GRID  if  it  is 
to  be  saved.  Valid  entries  arc  8,  9 or  10.  Default 
=0  fno  GRID  file) . 

F’lM.’H  GRID 

= 1 if  grid  file  is  to  be  piinchcd  (CIK!  filc=TAI’i;7) . 
= 0 otherwise. 

V.\RI.\BIdi  F-ORNIM  F'L\G  - for  Input  l>T>c  1.  Ill  is  is 
the  flag  to  indicate  whether  the  data  records  arc 
to  be  read  in  a particular  fonrnit  or  in  free  foniut . 
= 0 free  fomuit 
= 1 variable  format 

The  variable  format  is  input  on  the  following  card 
in  6A10  fonnat . It  must  be  enclosed  in  parentheses. 
•An  AZ  specification  must  be  included  as  eitlicr  the 
first  or  last  field.  Iliis  will  be  the  position  of 
tlic  two  nines  (99)  entered  in  the  last  record  of  the 
input  file. 


(PI.1ASI:  NOn,;  Using  the  Will.  GRID  or  UTtl  P.  GOl'.I-l'iniXrS  one  must 

remember  the  joli  number  of  the  step  as  it  is  printed  out  or  the  data 
c;innot  be  read  back  again.) 


Processing  I’arameter  Card 

Ihis  card  is  used  to  select  the  par;u:ieters  used  in 
determining  the  tlONIVHIR  GRID.  Ijitries  not 
specified  will  take  the  standard  svstem  default 
value. 

Panimcter  1 ; 4 

P.ir.'unetcr  Z:  NUMltl.R  (''1-  (X!I7V\T.S  - tlie  minimum  number  of  octant^- 

around  each  DATA  IXUNT  whicli  must  be  filk\i  with 
MDJIBOR'S  to  ensure  the  DMA  lAUNTS  inclusion  in 
the  contouring  process.  The  default  is  .5. 


r 


Parameter  3: 


Parameter  4: 


Parameter  5: 


Parameter  6; 


WPiarriNG  factor  for  TUF  data  points  - the  factor 

to  be  used  in  weighting  a NFIGIIBOR'S  contribution 
to  the  PLANE  surface  fitted  through  the  DATA  TOINT. 
The  default  value  of  -2  is  approximately  equivalent 
to  weighting  according  to  the  inverse  square  of  the 
distance  between  the  DATA  POINT  and  its  neighbor. 
NUMBER  OF  NEIGHBORS  - the  number  of  closest 
NEIGHBORS  to  a GRID  POINT  which  will  be  used  for 
determining  a value  at  that  GRID  POINT.  The 
default  is  10  and  the  maximum  permissible  is  20. 
WEIGHTING  FACTOR  FOR  GRID  - the  factor  to  be  used 
in  weighting  contributions  to  the  GRID  POINT  value 
from  the  PLANliS  of  the  neighboring  DATA  TOINTS. 

The  default  value  is  -2. 

RANGE  - the  radius  in  inches  of  the  search  area 
around  a GRID  POINT.  ITie  NIJMBF.R  OF  NI.IGHBOR's  must 
be  found  within  this  radius  if  the  GRID  POINT  is  to 
be  evaluated,  llie  default  is  5.0. 


Plot  Parameters  Card 

Users  may  alternate  3 pens  in  one  job.  Colors 
green,  red,  blue  and  black  arc  available,  llie 
basic  setup  for  the  Calcomp  93b  is: 

Pen  1 - black 
Pen  2 - blue 
Pen  3 - red 

A change  in  pen  arrangement  will  increase  the  basic 
charge  rate.  If  the  various  pen  options  are  used 
(other  than  1)  one  must  he  certain  to  give  tlie 
plotter  operator  the  appropriate  instructions. 
Parameter  1 : 5 

Parameter  2:  WTITIH  - the  rruiximum  width  of  paper  (in  inches)  to 

he  used  when  drawing  the  m<ap.  If  the  site  of  map 
is  greater  th;in  this  width,  the  map  will  be  drawti 
in  as  many  segments  as  it  rc()uires,  each  of  which 
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is  the  specified  width.  /\pproximatcly  1.5  inches 
is  required  for  title  information  and  so  for  11- 
inch  paper  a value  of  9 is  usually  optimal.  The 
default  value  is  28. 

Parameter  3:  PLOT  Ti  rLT; 

= 0 no  title  plotted  (default) 

= 1 through  3 if  title  is  to  be  plotted  with 
corresponding  pen. 

Parameter  4;  PLOT  .VXIS  - 0 to  5 have  same  significance  as  PIDT 
TITLL. 

Parameter  5:  PI.OT  SiTIBOLS  - 0 to  3 as  above. 

Parameter  b:  Si'MBOI.  NlJNffiLR  - any  noncero  number  will  select  the 

corresponding  plotter  s\mibol  (see  SYMBOL  TARLl;  in 
the  plotter  manual)  to  be  drawii  at  the  X,V  location 
of  each  data  point.  Itcfault  is  the  svinbol  X. 
Panuncter  7:  .S'l’MBOL  Sill.  - height  of  SNTnbols  in  inches,  default 

is  .o:". 

Panimeter  8:  PLOT  DATA  - 0-3  with  s;imc  significance  as  I’LOT 

TITLL.  Default  = 0. 

Parameter  9:  DATA  SIZL  - the  height  in  inches  at  witich  the  data 

values  will  be  drawn  if  PliYl'  DATA^O  (Panimeter  S). 
Default  = .0"' . 

I’animcter  10:  DATA  lYlSITlON 

=0  if  posted  data  values  are  to  be  centered  on  the 
DATA  IXl ! NT  position. 

=1  if  values  to  be  posted  to  the  right  of  the  DATA 
(XlIM  . 

= 2 for  above  the  point. 

= 3 to  the  left  of  the  point. 

= 4 for  below  the  data  point. 

Parameter  11:  DATA  .WC.I.I:  - tlic  angle  in  degrees  at  wliich  data 

values  are  to  be  posted.  Hiis  ;ingle  is  measured  in 
degrees  from  the  positive  \-a\is.  l)t'fault=0. 
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Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Parameter 

Par;imcter 


12;  DATA  FOI^T  - the  number  of  decimal  places  to  the 
right  of  the  decimal  point.  Default=0  (none). 

13:  PLOT  HEAVY  CONTOURS  - pen  number  1-3  or  0 for  no 

plot.  Same  as  PLCTT  TITLE.  Heavy  contours  are 
plotted  twice  to  give  a bold  line. 

14:  PLOT  LIGHT  CONTOURS  - 0-3  as  above. 

15;  ANNOTATE  HEAVY  CONTaiRS  - a value  0 to  3 as  above. 

If  nonzero,  the  numerical  value  of  the  heavy' 
contour  lines  will  be  plotted  along  the  contour  line. 

16:  ANNOTATE  LIGHT  CONTOURS  - same  as  above  (Parameter 

15)  except  for  light  contours. 

17:  AN7MOTATION  SIZE!  - the  height  (in  inches)  at  which 

contour  annotations  are  to  be  drawn.  Default 
value  =0.1  inches. 

18;  ANNCT’ATION  FORM\T  - number  of  decimal  places  to  the 
right  of  the  decimal  point  for  contour  annotation. 
Uefault=0  (none). 

19;  HEAVT  CONTOUR  FREQUENCY  - a number  \ such  that 
every  Nth  contour  will  be  drawn  heavy. 

20;  MINIMUM  AN'SCTiATlON  LENITTH  - a value  in  inches  such 
that  any  contours  shorter  than  this  value  will  not 
be  .'innotated.  l)cfault=5  inches. 

21:  MINIMUM  NUMBER  OF  POINTS  - a value  sucli  tlvat  any 

contour  based  on  fewer  defining  points  will  not  be 
drawn,  llie  default  value  is  4.  For  maps  to  be 
drawn  in  more  than  one  segment  this  par;imcter  should 
be  specified  cc|ual  to  1. 

22;  MINIMUM  I.ENimi  OF'  CONTOURS  - a value  in  indies  such 
that  any  contours  which  arc  shorter  than  this  value 
will  not  be  drawn.  Hie  default  value  is  3 times  tlie 
GRID  SIZE  specified.  If  maps  are  to  be  drawn  in 
more  tlian  one  segment  tliis  parameter  should  be  set 
to  0.01  in  order  to  ensure  matching  of  contour 
lines  at  the  boundaries. 


Parameter  23; 


CONTOUR  LOWliR  LIMIT  - any  contours  below  this  value 
will  be  omitted. 

Parameter  24;  CONTOUR  UPPLR  LIMIT  - any  contours  above  this  value 
will  be  omitted.  If  both  contour  limits  are  zero, 
the  contour  limits  are  chosen  from  the  data. 

Parameter  25;  MM*  SLPjXRVITON  - the  contour  map  will  be  plotted  in 
the  normal  fashion  up  to  and  including  this  contour 
level.  Tlien  a new  map  will  be  started  beginning 
with  this  contour  level,  lliis  is  the  feature 
called  M\P  SLP(\RVnON  discussed  in  Section  II.  If 
the  following  two  parameters  arc  ..c-o,  all  contours 
are  drawn  on  one  m;ip. 

Parameter  2b:  PLOT  IllcWT  CONTOURS  - same  as  Par;unetcr  13  except 

applies  to  separated  map. 

Parameter  2";  PLOT  LKiri  CONTiXIRS  - as  Parameter  14  for  separated 
map. 

Parameter  28;  .WNOiAlT  lIlcUT  - as  Parameter  15  for  separated  map. 

Parameter  29:  .VN'NOTATl . l.hlin  - as  Parameter  lb  for  separated  map. 


lixtra  Options  C3ard 
i’arameter  1 : 
Parameter  2: 


Panimctcr  3: 


(1 

SUPPRIMl'OSrnON  PLAG  - flag  to  define  a super- 
imiiosition  of  one  map  on  ;mother. 

- 0 no  superimposition 
= 1 superim|iose  tins  map  with  the  next. 

Default  value  is  0. 

NOISi;  TOl.I.llVNCL  - a value  which  determines  when  a 
contour  has  no  real  meaning.  Mien  two  Z-values  on 
the  grid  are  less  th;ui  or  equal  to  this  value  in 
distance,  tliat  contour  will  lie  eliminated.  I\?fault 
value  is  0.0001. 

RVSII  LLNCni  - real  number  to  i.lefine  the  length  in 
inches  of  a dash  or  space  in  the  dashed  lines. 
Default  value  is  0.0b. 


Parameter  4: 


Parameter  5: 


D/VSH  I’A'ITliRN  - integer  number  to  define  the  pattern 
of  the  dashed  line  contours.  Default  value  is  0 (no 
dashed  lines).  ITie  dash  pattern  is  a 12-bit  pattern 
with  each  bit  representing  a one  dash  length 
interval  along  the  contour. 


Dash  Pattern  Value 
2730 

3510 

3789 

3900 

4030 

4094 

4092 

4032 


Dash  Pattern  Value 
4078 

40-4 

4010 

3818 

3902 

3-54 

32~0 

2340 

38-0 


— 4044  _ 3970 

Parameter  O:  JUNIT  - tape  number  where  outline  points  are 

located.  (Xitlinc  points  may  define  some  body  or 
detail  to  be  superimposed  on  the  contour  map. 
Default  value  is  0 (for  no  outline  points). 

= 5 for  card  input 

= 8,9,10  for  pcniianent  file  or  tape  input 


A card  with  HNDMAP  tv^icd  in  AO  format  must  appear  Ix'twoen  e.Kh  t'.' 
sets  of data  cards  to  produce  a contour  map. 


r.nd  of  .lob  Tire  last  l.NDM\l’  card  of  the  last  set  of  data  cards  mav  be 

followed  by  a card  witli  a ”9”  in  col.  1,  or  by  a -/8/9  "end  of 
record"  card  (or  0/7/8/9  "end  of  file"  card  if  the  input  file  is 
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2H 
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other  than  the  card  readerl . 


INRn-  DATA  CARDS 

Following  is  an  explanation  of  the  preparation  of  the  five  different 
input  types. 

Input  Type  1 - X,Y,Z  Coordinates 

This  data  may  be  input  in  free  format  or  a particular  format 
may  be  sjiecified  by  setting  the  flag  in  Parameter  lb  of  the 
Input/Output  card.  In  either  case,  data  consists  of  one  set  of 
X,Y,Z  coordinates  per  record.  The  series  of  input  data  prints 
must  be  terminated  with  a record  containing  the  characters  99. 
,\ny  X,Y,Z  values  appearing  on  the  99  card  are  ignored.  Special 
care  must  be  tiiken  if  the  free  format  is  used.  Free  format  will 
always  expect  to  read  four  parameters  per  record  unless  a 
slash  is  encountered.  Tliereforc,  each  data  record  must  read 
X,Y,Z,0  or  X.Y,Z/  and  the  end  record  must  read  0.,0.,0.,99;  or, 
depending  on  the  end  test  position  indicator,  0,\,Y,Z  for  each 
data  record  ;ind  99,0.,0.,0.  or  99/  for  the  end  record.  It  imiv 
facilitate  the  free  format  input  to  indicate  that  the  99  will 
he  the  last  characters  of  the  record  (i.e.,  I’arameter  9=1  on 
Input/(\itput  card). 

Input  Type  2 - Trend  Surface  Coefficients 

ITiis  t>’pe  of  data  is  read  in  free  format.  Hie  data  record 
consists  of  the  order  of  the  trend  surface  (1,2, 3, 4),  the 
number  of  coefficients  (3,b,10,15,  respectively)  followed  by 
tlic  coefficients.  No  end  of  data  specification  is  required. 

Input  r>T>e  3 - (irid  Points 

lliis  tvj)c  of  data  is  read  in  free  fonnat . The  first  record 
contains  the  number  of  rows  and  the  mmiber  of  columns  in  the 
grid.  F.ach  succeeding  record  (as  rruiny  as  required)  contains  up 
to  seven  grid  point  values 


Input  4 vid  5 


These  are  formats  internal  to  the  program  produced  using  the 
appropriate  read  and  write  options.  Input  option  TYPI->5  must 
have  parameter  5 of  card  5 equal  to  0 (plot  symbols)  and 
parameter  8 of  card  5 equal  to  0 (plot  data)  or  the  program  will 
fail . 

OUTLINE  INPUT  DATA  CARDS 

Outline  points  are  specified  as  X,Y  coordinate  pairs.  They  arc  read 
in  free  format,  one  pair  per  record.  To  lift  the  pen  (specify  a "hole") 
the  user  inputs  998.,  998.  . The  end  of  data  indicator  is  999.,  999.  . 


.■^0 


V.  SAMPLE  RUN 


Below  is  a sample  of  a simple,  straightforward  CONTOUR  run.  One  map 
was  produced  from  one  set  of  X,Y,Z  coordinates.  The  input  was  as  follows: 

CAEMCT ,016000 ,T180 ,P3 ,MT1 . 1843/MARQUARDT 

Q L\RGE , ( j\EM , nnnnnnnnnn . 

ATTACH , CONT , CAMVCONTOURTSK , I D=GiMV . 

VSN(TAPE11=CB0536) 

REQUEST, T/\PE11 , HI  ,RING. 

CONT. 

789 

1 Tl-ST  CONTaiR  M\P 

2.100. . 5. .100. .5. .5.  ,21, .4, 21  , .4 

3. 1.5. 0. 0. 0.0. 0.1. 1.1. 1.0. 0.0.1 
(F3.0,2E4.0,A21 

4. 3.0, :’, 0,10. 

5.8. . 1.1. 1.2.. 02.1.. 07. 0.0.  ,0 , 1 ,1  ,1  , 1 , . 1 ,0 ,5 , 10*0 
6,5*0 

148  125  140 
122  131  145 
143  145  152 
158  149  146 
115  155  155 
135  15~  i:’3 
155  165  156 
185  168  145 
165  170  15:' 

135  172  154 
112  181  145 
^5  182  170 
153  183  143 

149  186  140 
182  193  145 

99 

i;.\nMVP 

9 


Iliis  job  consists  of  15  data  points  having  the  following  mininvi  ;uid 
maxima. 

Minimum  MxLximum 
X 112  185 

Y 125  193 

: 140  r3 
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Both  the  X origin  and  Y origin  have  been  chosen  as  equal  to  100.  The  ma.x - 
imum  required  for  both  X and  Y is  200.  A grid  size  of  5 data  units 
corresponding  to  .4  inches  has  been  selected.  The  range  of  both  X and  V 
is  100  data  units,  therefore,  at  5 units  per  grid  square,  20  squares, 

(X  and  Y in  this  case)  are  required.  Hence,  the  number  of  columns  and 
number  of  rows  in  the  grid  matrix  both  equal  21  for  a total  of  441  ('.RID 
points  (21  X 21).  The  map  size  is  .4  x 20  (both  directions)  or  8 inches 
by  8 inches.  The  input  is  from  card  images,  hence  the  input  unit  is  5. 

The  order  of  the  coordinates  is  X,Y,Z  and  the  input  format  is 
(F5.0,2I-4.0,A2) . The  A2  is  for  the  99  value  in  the  last  record  used  to 
signal  the  end  of  the  data.  Output  requests  are: 
a list  of  the  input  data. 

a list  of  the  coefficients  of  the  PLANTS  generated, 
a list  of  the  GRID  matrix  created. 


.•S2 


JOB 


1 


i 


TEST  CONTOUR  mkP 
«iP 

SIZE  IN  0PIGIH4L  04T4  JNI»> 


X - 

100.080  UNITS 

TO 

ZOO. 000  UNIT^ 

G R 1 1 S T ,■»  - 

000 

UN  I T s 

V - 

100. on  UNITS 

TO 

zoo. 000  UNITS 

GRIT  SIZE  - 

S.  000 

units 

pmtsical 

MAP 

SIZE 

X - 

A. 

00  INCHES 

GRID 

SIZE  .<.•  INCH 

▼ - 

8. 

00  INCHES 

^R  13 

SIZP  .-0  INCH 

SIZE  IN  CONPUTEB 

Zl  columns  9r  Zl  pokS  TOm  **41  SPIO  SQUktfES 


INPUT  iNOlZMTtO  - *,T,2  Okll 

input  unit  5 

REWtNO  0 
INPUT  JOB  0 
X VftRIABLE  1 
T VARIABLE  2 
Z VARIABLE  3 
pos.  or  eno  test  i 


!JT  PftP*N£T£R5 

IN  SCO 

PRINT  -^AT*  I 

print  CI'PPICHNTS  I 

PRINT  GRID  I 

OUTPUT  COer  ON  UNIT  Q 

OUTPUT  GRID  ON  UNIT  0 

OUTPUT  GRID  ON  CAROS  0 

INPUT  PORNAT  (P3.0»0,A2) 


HUMBER  Of  OCI»NTS  NMICH  >1UST  3E  FILLET  E F »L  J»  T I •J'l  0»T»  POINT  PL»NES 


NdjlKd  hS«09I*  i NiOM31  HSfO 


exccurioN 


HESS*&fS 


REI01H&  I,t*2  0»TA  XN  BCD  - START  . ,223  SECONTS 

END  * . 2?B  SECONDS 

- TOTAL  * .aur  SECOHOS 

NO.  Of  DAT*  POINTS  » IS 


} 


COORCINATES  Of  CONTOURING  DAT* 


X 

T 

2 

149. OOOOOODO 

125. OOOOOOOD 

140. 00001000 

122,  00000090 

1 11. 00 JOOOOO 

145. ooaooooo 

1W3. OOOOOOOO 

1W5. OUOOOOOO 

152. OOOOOOOO 

159. 00000000 

l->9.  OOOOOOOO 

146. OOOOOOOO 

U5.00Q09Q91 

155. CQ900C0Q 

155. aooooooo 

135.00000000 

157. 00000000 

173. OOOOOOOO 

15S. 00000000 

165. OOOOOOOO 

156. OOOOOOOO 

165.0C000990 

169. C0300C00 

145. OOOOOOOO 

165. 00000000 

1 70. OOOOOOOO 

157. OOOOOOOO 

135.00000000 

172. OOOOOOOO 

154.00000000 

1 1?. OOOOC090 

191.00000000 

1**5.  OOOOOOOO 

175. 00060000 

192. OOOOOOOO 

170. OQOQOOOO 

153.00000000 

191. OOOOOOOO 

141. OOOOOOOO 

149. OOOOOOUO 

196. OOOOOOOO 

140. OOOOOOOO 

102. 00000000 

191.  OOOOOOOO 

145. nooooooo 

E*ECO^IOn  NESSAGfS 

Evaluation  of  planes  through  oata  points  - st4r» 

s J 

- total 


.S?’  SF^-'NT' 
. T Sf 
.0?»  .ECPNCS 


COEFFICIENTS 


P L » H E S 


Through 


C(l) 

CC2I** 

C(3»*T 

. 370?6{.TE*03 

•13333336*00 

-.21000106*11 

. W8194l2r»l3 

.1176471£»00 

.62941186*00 

• 1701725E»I3 

-.07123496*00 

• 73389006*00 

• lt»2«*9**Sf*0  3 

-.74242046*  00 

.67657116*00 

8961S3«C»03 

. 15304626*0 1 

• 30769236*  00 

• 1670?S7E»03 

-.54660106-01 

• 6 5 2 0 6 4 86  - 01 

• 16403S8E»03 

-.61995636-01 

.95362786-02 

-.292$57ie*0l 

-.75714296*00 

. 17142866*01 

• 73l67a3’r»02 

.^2260466*  00 

-.20763316*00 

•363696ie»03 

-.20036006*00 

-.10619046*01 

• 1532613E»03 

-.23437«0E-0l 

-• 31250006-01 

•24€990aE»03 

-.22006006*00 

-.20917736*00 

. 22664056*03 

.46730026*00 

-. 84701266*00 

.29909426*03 

-.13096236*00 

-.74«.026D6*00 

-.29560006*03 

.30000006*00 

•20000006*01 

Execution 


N E s s « G E s 


DETERNINRTIOM  of  contour  GRIO  - ST»RT  I .k?i,  SrcONOS 

end  » SECONDS 

- TOTXL  = l.bSS  SECONDS 

NUMBER  OF  grid  points  EVALUATED  = GNl 
CONTOURING  LIMITS  CHOSEN  - LONER  = .lN000E*e3 

- UPPER  = .lT!00f«03 
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IMUJUIUJiyWWUiUJteiWMWWWMMlWMWWMMMUflWWWUMMMlWMMMWWWWiyWlMlil 
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lAi  tM  W tei  Ui  M W 1^  ^ ^ 1^  ^ U M VJ  W W W M ^ yj  1^  ^ 14/  |^^|^g^^|4/|^|4/|4/4^^^l^^|^ 


^ ^J  ^ yj  ^ ^j  ^ yj  ^ ^ ^ ^ ^ ^ ^ 


o 


1 

f 


o 

(k 


ItiLt.  LAjtxJLA.L«Jt^jL4/LtJL..  L.^y»LAj'j^itAjLAjL»JtAjL*.’L^UJU-«tA^MJU-.»Aji^t»-UJt^t^i^L>-tAjiAjUJL^L^IJ^L«^L».t^^l^ 

oc«ou^Tir«  #(f'^tT'aw'f^c‘^iro'-^<^4r«A(*«A  .#9a<^<^v«u'^X^ 


oeoeo9ooeoooe'o‘9r99ec99eeaooa0oooeaaoa9«3eooaeoe 


yj  ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ IM  ^ ^ ly  ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ W W W ^ ^ ^ ^ 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

(1)  DTNSRDC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

(2)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
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